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THEORETICAL SYMMETRIC SPAN LOADING DUE TO FLAP DEFLECTION FOR WINGS ОҒ 
ARBITRARY PLAN FORM AT SUBSONIC SPEEDS! 


Ву Јонх DEYounG 


SUMMARY 


A simplified lifting-surface theory is applied to the problem of 
evaluating span loading due to flap deflection for arbitrary wing 
plan forms. With the resulting procedure, the effects of flap 
deflection on the span loading and associated aerodynamic 
characteristics can be easily computed for any wing which 18 
symmetrical about the root chord and which has a straight quar- 
ter-chord line over the wing semispan. The effects of compressi- 
bility and spanwise variation of section lift-curve slope are 
taken into account by the procedure. 

For the case of straight-tapered wings, the load distribution 
and lift due to flap deflection are presented in chart form for a 
broad range of wing plan forms. Experimental values of flap 
effectiveness are compared with the values obtained from these 
charts for several wings. The comparison shows good agree- 
ment between theory and experiment. 

The method presented can also бе used to calculate the down- 
wash in the vertical center of the wake of a wing which has 
arbitrary spanwise loading. 


INTRODUCTION 


Reference 1 presents a procedure which is based on simpli- 
fied lifting-surface theory and which permits the rapid 
determination of the symmetrie spanwise distribution of lift 
on wings having essentially arbitrary plan form. With 
reference 1, the symmetric spanwise loading can be found for 
any symmetrie angle-of-attack distribution such as that due 
to spanwise continuous twist or camber. Тһе method satis- 
fies the boundary condition at only a finite number of span- 
wise points; thus, to insure accuracy for the loading obtained 
by this method, the angle-of-attack distribution must be 
continuous and must not vary excessively. If the angle-of- 
attack distribution is discontinuous, such as that due to 
partial wing-span flaps, or varies excessively, the boundary 
condition should be satisfied at sufficient points to insure а 
correct answer. This normally entails an excessive amount 
of labor, since the computations required increase approxi- 
mately as the square of the number of points used. 

An alternative procedure is to use a method which satisfies 
the boundary condition at only & few points where this 
boundary condition is adjusted so that at these points the 
calculated loading is identical with that which would be 


1 Supersedes NACA ТМ 2278, “Theoretical Symmetric Span Loading Due to Flap Deflec- 
tion for Wings of Arbitrary Plan Form at Subsonic Speeds” by John DeYoung, 1951. 








predicted by a method satisfying the boundary condition 
at all points. This specific boundary condition for the case 
of antisymmetrically deflected ailerons is developed in refer- 
ence 2, which is based on the zero-aspect-ratio theory of 
reference 3. А similar development can determine the 
specific boundary condition for the case of symmetrically 
deflected flaps. 

It is the purpose of the present report to determine the 
specific boundary condition and to provide simple methods 
for computing the aerodynamic characteristics for wings with 
symmetrically deflected flaps. 


NOTATION 
. (5? 
а aspect ratio (5) 
Aya coefficient depending on wing geometry and 


indicating the influence of symmetric load- 
ing at span station n on the downwash 
angle at span station v 


b span of the wing measured perpendicular to 
the plane of symmetry, feet 

с wing chord,? feet 

Cas mean wing chord? (2) feet 

€ local lift coefficient (к=нг 

Cr lift coefficient ca) 

С, rate of change of lift coefficient with flap 


deflection, per radian 
Ce rate of change of lift coefficient with wing 
angle of attack, per radian 





ыр spanwise loading coeffieient for unit lift 
C'rlar 
š (226) 
coefficient | -一 一 
Cr 
а, scale factor 
ба: factors of loading interpolation funetion 
G spanwise loading coefficient or dimensionless 
5 a Cic E 
" circulation ($ Э ог (; ғ) 
6 spanwise loading coefficient per radian of 
flap deflection (È) 
1 


2 Measured parallel to the plane о? symmetry. 
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wing-chord distribution parameter (at span 
station >) which includes the effects of 
compressibility and section lift-curve slope 


[C] 


integration factors for spanwise loading due 
to flaps 

Mach number 

arbitrary number of span stations 

free-stream dynamic pressure, pounds per 
square foot 

wing area, square feet 

free-stream velocity, feet per second 

induced velocity, normal to thelifting surface, 
positive for downwash, feet per second 

longitudinal coordinate measured from the 
lateral plane through the quarter chord of 
the wing-root chord ? feet 

lateral coordinate measured from the wing 
root perpendicular to the plane of sym- 
metry, feet 

wing angle of attack? radians 

section angle of attack at span station 23 
radians 


compressibility parameter 4|1 — M?| 

circulation, feet squared per second 

angle of deflection of Пар, radians 

angle of deflection of flap measured perpen- 
dicular to the hinge line, radians 





Сі 
lift-effectiveness parameter? (5 š 
s La 


: ; : “di y 
dimensionless lateral coordinate (she) 


dimensionless fap span on one wing panel, 
measured perpendicular to the plane of 
symmetry, from the wing root outboard 
for inboard flaps, and from the wing tip in- 
board for outboard flaps (mm = 

trigonometric. spanwise coordinate defining 
the spanwise end of the flap, radians 

ratio of section lift-curve slope at span sta- 


tion ў to = both at the same Mach 


number 4 

sweep angle of the wing quarter-chord line, 
positive for sweepback, degrees 

sweep parameter which includes the effects of 


gita {ап А 
compressibility (tan =) degrees 


wing taper ratio (Sp chor’ 
SSR © \root chord 


dimensionless longitudinal coordinate (52) 


trigonometric spanwise coordinate (cos! n), 
radians 


* See footnote 2, page 1. 
* See footnote 2, page 1. 

















SUBSCRIPTS 
n integer pertaiuing to span station associated 
. | пт 
ith the loadir (- =) 
with the loading ( ў —cos EET 
y integer pertaining to span station associated 
with the downwash (п--соз LR т) 
k half-integer pertaining to span station 
( =cos т 
ща т--1 
f pertaining to Нарв 
i wing tip, 
ғ wing root 
ар average ог mean 
1 denoting full wing-chord Варз (“= 1) 


DEVELOPMENT OF METHOD 


The simplified lifting-surface theory used herein replaces a 
lifting surface by a lifting vortex located at the wing quarter- 
chord line. The boundary condition specified for determin- 
ing the vortex strength distribution is that, along the three- 
quarter-chord line of the wing, there shall be no flow through 
the lifting surface. 

The mathematical method of solution and discussion of the 
limitations of the theory are given in detail in NACA Rep. 
921, 1948 (reference 1). The same mathematical procedures 
and assumptions are used in the present report. 

The mathematical solution for the determination of span- 
wise loading consists of satisfying the boundary condition 
that along the three-quarter-chord line of the wing there shall 
be no flow through the lifting surface. If the boundary 
condition (i. e., the spanwise twist along the three-quarter- 
ehord line) can be approximated by a few terms of a mathe- 
matical series, then the “no-flow” condition need bo satisfied 
exactly at only a few spanwise stations and the calculated 
loading can be expressed by a few terms of a mathematical 
series. Thus it has been found that satisfying the boundary 
conditions at four spanwise points enables accurate compu- 
tation of the loading due to symmetric twist (reference 1), 
whereas only three spanwise points are required for accurate 
computation of the loading due to antisymmetric twist 
(reference 2). Тһе distribution of spanwise twists in these 
two cases is such that at most a four and а three term trig- 
onometric series, respectively, are needed to duplicate Шеш. 

For the case of irregular and discontinuous distribution of 
twist such as that required to simulate the span loading due 
to deflection of a partial-span flap, it becomes necessary to 
satisfy the boundary condition at a large number of points, 
particularly near the point of discontinuity, in order that 
sufficient accuracy will be insured. Since the computation 
labor inereases approximately as the square of the number of 
points considered, the foregoing method becomes less ргас- 
tieal. An alternate approach is to alter the boundary con- 
dition for the method involving few points at which the 
boundary condition is satisfied, so that at these points the 
loading is identical with that obtained by the method for 
which the boundary condition is satisfied at all points. 
This process is practical because the spanwise loading of 
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wings becomes independent of plan form as aspect ratio 
approaches zero (reference 3). Thus the solutions for the 
zero-aspect-ratio wing which satisfy the boundary condition 
at all points can be used to obtain the proper adjustment of 
the boundary condition for a method that satisfies the bound- 
ary condition at only a few points. With the altered bound- 
ary condition, the method satisfving conditions at only a 
few points will give accurate results for plan forms of very 
low aspect ratio, and should give good accuracy for wings 
of high aspect ratio since the alteration of the boundary 
condition can be kept small by proper choice of flap spans, 
and also since loading distribution is affeeted less by plan 
form than by the boundary condition. The mathematical 
procedure for determining the altered boundary conditions 
is developed in appendix A. 

Although an altered boundary condition can be found for 
an arbitrarv position of the point of discontinuity, the 
alteration can be made a minimum and of small order by 
locating the discontinuity at the mean position between 
two points at which the boundary condition is satisfied. 
With the point of discontinuity defined, the flap span be- 
comes defined. The results for arbitrary Нар spans аге 
found by cross plots of results for the defined spans. 

When the flap spans and corresponding twists (altered 
boundary conditions) have been determined, the method is 
applied in the same manner as that of reference 1. In this 
way, (m4-1)/2 linear equations in terms of symmetric- 
loading distribution are obtained whieh satisfy the wing 
angle-of-attaek conditions ° at the three-quarter-chord line 
at m stations, where m is ап arbitrary odd integer. Тһе 


application in the present report is with m=7. The linear 
equations are represented by the summations 
4 
y= >) din Си, p—1,2,8,4 (1) 
nat 


where the integer subseripts у and n pertain to span stations 
n=cos (рт/8) and cos (пх/8), respectively. To obtain the 
loading coefficients G, from the simultaneous solution of 
equation (1), it remains only to evaluate the coefficients 
п, (which depend on wing geometry, compressibility, and 
section lift-curve slope) and the spanwise variation of the 
symmetric angle of attack а,. The reader should note that 
with а, constant (spanwise) equation (1) represents additional 
loading, and that with а, varying symmetrically and con- 
tinuously it represents basic loading (provided the net lift 
is zero). 
EVALUATION OF ара 


As shown in reference 1, а,, can be conveniently presented 
аз a function of two parameters,® namely, a sweep parameter 
defined as Ag=tan7! [(tan А)/8 and a wing-chord distribution 
parameter H, defined by 


H,=d, (=) (=a) 


5 As discussed in reference 1, e, is not limited to small values but can be as large an angle 
as desired, provided separation does not occur. 

š The effects of compressibility and section lift-curve slope are equivalent to a change in 
wing plan forms and can be taken into account by a proper adjustment of the а», values. 





(2) 











where 
к, ratio of experimental section lift-curve slope 
at span station v to the theoretical value of 
= both at the same Mach number 
с, wing chord at span station > 
d, scale factor which has the following values: 
0.061 for »=1 
234 for y —2 
.381 for »—3 
320 for у=4 


Equation (2) can be written in an alternate form that gives 
H, in terms of wing geometry parameters that are more 
significant; thus 





‘BA 1 
Didi (5 ) я 3 
Каз (куј кар) (¢y/Caz) ( ) 
where 
К ratio of average section lift-curve slope to 
27/8, both at the same Mach number 
ua spanwise distribution of section lift-curve 
^et slope for a given Mach number 
с, ORE а лова Ди : 
= spanwise distribution of the wing chord 
‘ay 
A : 
ва aspect-ratio parameter 
Е 


The term 1/(к/кол) (с Са) of equation (3) gives the effective 
aerodynamic taper of a wing. The distribution of к/к, may 
vary with Mach number, particularly at transonie speeds 
(e. g., due to spanwise variation of airfoil section). However, 
since this distribution contributes to taper effect, the primary 
change will be in the loading distribution rather than in gross 
load. 

With H, determined from equation (2) or (3), the values of 
аул, 16 in all, can be found in figure | where а,, is given asa 
function of Я, for various values of Ag. 

For the case of straight-tapered wings with arbitrary sec- 
tion lift-curve-slope distribution for which the chord distri- 
bution is specified by taper ratio, evaluation of equation (3) 
is facilitated through the use of figure 2, which gives к, ,/8.А 
for each of the four span stations as a function of taper ratio. 


EVALUATION OF as 


The symmetric angle-of-attack distributions most com. 
monly encountered are those resulting from constant angle of 
attack, twist and camber, and the effective angle of attack 
due to Нар deflection. The first two cases are the subjects 
of reference 1. The effective angle-of-attack distributions 
due to flap deflection are the subject of the present report. 

To minimize the computations involved, it is convenient 
to consider the cases of inboard and outboard flaps inde- 
pendently. The inboard flaps are measured from the wing 
midspan outboard, and the outboard flaps are measured 
from the wing tip inboard. Three flap spans for inboard 
and outboard flaps are considered in appendix A and the 
equivalent angle-of-attack distribution for unit flap deflection 
a,/6, for each case is listed in table A4. The values of a,/6, 
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given in table А4 are for (ej/c) 21. The results can be easily 
modified for the case of (с/с)<1 as wil be shown later. 
Тһе ease of full wing-span flaps is easily considered since 
this is the ease of constant angle of attack which results in 
additional-type loading. 


APPLICATION OF METHOD 


For the cases of symmetric angle-of-attack distributions 
over & wing with flaps deflected, itis possible to form the set 
of simultaneous equations which are required for the solution 
of loading distribution for an arbitrary plan form. With the 
loading known, integration formulas can be given to deter- 
mine aerodynamie coefficients. 

'The loading distribution G,, determined from the solutions 
of the simultaneous equations, is a function of a,,, which has 
been shown in a preceding section to be а function of wing 
geometry, compressibility, and section lift-curve slope. The 
aerodynamie coefficients are obtained by integrations of the 
laoding distribution and, therefore, will also be functions of 
wing geometry, compressibility, and section lift-eurve-slope 
parameters. 

In the following sections, application of the method to 
obtain the general solution for arbitrary spanwise distribu- 
tion of wing chord is outlined and solutions are presented for 
the special case of wings having straight taper. 


GENERAL SOLUTION 


Simultaneous equations for the logding.— The а), values 
to be substituted in the equations that follow are obtained 
from figure 1 and table 17 with values of H, given by equation 
(2) or (3). 

1. Deflected inboard flaps: The flap spans measured from 
the wing midspan outboard are given by т. The set of 
simultaneous equations for symmetric spanwise loading due 
to deflection of any of the three following flap spans can be 
obtained from the appropriate set of the following equations: 





























Case I H | ш 
| 1 » 
W 0.195 | 0.556 | 0.831 
= —0.017 | 0.032 м =0 оо азбан 
4 09 | —.003 | 1.021 | =алб фор јаз Рамба 
й | е жу на ы AC) 
ш —.014 .994 958 | =ап азба) ан 
si . 988 . 976 1.010 хабағаабығаабы-анб; 
я : 
where 


G,= 6.8. 


For the case of ту=1; a,/è; is unity over the span of the wing 
(same condition as solution for additional loading). 

9, Deflected outboard faps: The Вар spans measured from 
the wing tip inboard are given Бул, The set of simultaneous 
equations for symmetric loading due to deflection of any of 
the three following flap spans can be obtained from the ap- 
propriate set of the following equations: 


? Values of ae, beyond the scope of figure 1 and table I can be obtained from equation (B2), 
which gives the linear asymptotes of the а, function. 























— n . 

Case IV Y VI 

я, 0.169 0.444 0.805 

e = = = = 

hr 1.011 0.968 1.017 | ап Фан ов Фоа 

= —.021 1.003 972 апан азбан, 

оз = = = = (5) 

= 045 008 1.014 жамбі-анба--ааба-ғазбі 

x | —.010 .024 .012 ха фаб атан, 
where 

G,= Gaf: 


For the case of 9-1, о,/$у=1. 

Loading distribution.— The spanwise loading distributions 
сап be determined for various flap configurations as follows: 

1. Full wing-chord flaps [(сује) 21]: The loading is zero at 
the wing tip and is determined at four other span stations of 
n=cos (пт/8) =0.924, 0.707, 0.383, and 0 by solution of the 
appropriate set of simultaneous equations from the previous 
section. With equation (А9) and tables А6, A7, and А8, 
the loading can be interpolated at span stations 7=0.981, 
0.831, 0.556, and 0.195 for each of the flap spans considered. 
With these given points and the knowledge that the slope 
of the loading-distribution curve is theoretically infinite 24 
the point of angle-of-attack discontinuity or flap spanwise 
end, the loading distribution can be faired. For flap spans 
other than those given, the loading distribution of flaps that 
are neither inboard nor outboard can be obtained from the 
concept that loadings are additive (discussed later). 

2. Constant fraction of wing-chord flaps [(e Jay = constant]: 
The spanwise loading distribution is equal to the product of 
the loading distribution due to full wing-chord flaps and the 
lift-effectiveness parameter 00/48. The factor da/di is a 
function of the ratio of flap chord to wing chord с//е. Values 
of da/dò are presented in figure 3, which is reproduced from 
figure 18 of reference 4+- Although reference 4 from which 
figure 3 was taken limits the Mach number range {о values 
less than 0.2, this limitation is believed to be too restrictive 
since linearized theory indicates that de/dé is unaffected by 
compressibility for the two-dimensional wing. As indicated 
in reference 3, da/dé is strongly affected by aspect ratio if 
the aspect ratio (ог 8.4) becomes very small; however, com- 
parisons with experiment indicate the values of de/dé from 
figure 3 appear to be valid for BA>2. This limit to which 
da/dé can be used will be considered further in the section 
Discussion when comparisons with experimental results are 
made. _ 225 = Ж. 

8. Arbitary spanwise distribution of flap chord (с/с)- 
variable]: The flap can be divided into several parts cach 
having constant 00/28 and the load distribution due to each 
part determined. The total distribution is then the sum of 
these individual load distributions. 

Lift coefficient.—The lift coefficients due to flap deflection 
for the same flap configurations previously discussed can be 
found as follows: 

1. Full wing-chord flaps [(c,/e)=1]: The spanwise loading 
due to flap deflection is, in general, so complicated that 
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numerical integration based on four spanwise values cannot 
be performed accurately with conventional integration 
formulas. However, in appendix A a special integration 
formula (which needs but four spanwise values) is developed 
which applies to the different flap spans. Equation (Аб) 
can be written as : 


8C 3 = = = = 
=== (85) MO 1 0-0, thee (6) 


where, for each of the cases of equations (4) and (5), Ше А, 
values are given by 








Inbosrd flaps 








Case 











лу 











== 
м 0. 301 0.301 0.301 
ha 556 555 
hi 128 .725 726 
ы 395 .393 393 











2. Constant fraction of wing-chord Faps [(сујс) =сопвіапі]: 
The flap effectiveness is given by 








(7) 


3. Arbitrary spanwise distribution of flap chord (сус) = | 


rariable}: Flaps for which ¢,/e varies spanwise on the wing 
can be considered as equivalent to a wing-twist distribution. 
The effective symmetric twist of the wing is given by 


da 
aL Š (8) 


where da/dé is now a function of spanwise position. If 
equation (8) is a continuous distribution such that it can be 
plotted by specifying its value at four semispan points, then 
the wing can be considered to be twisted and solutions found 
directly as for basic loading. 

When equation (8) is discontinuous, the angle of attack 
can be divided into spanwise steps of constant angle of 
attack and the total lift can be found by the summation of 
the lift due to each spanwise step. The lift of a spanwise 

ве 


L 
51 





step is obtained from a curve of lift coefficienti as a 


funetion of the extent of flap span from the wing root out- 
board. Such a curve can be obtained from the solutions of 
equation (6). Аз previously indicated for ту:=1 


(а ) Е ВС, 
Ка: pl Каз 


The lift coefficient due to Ше twist given by equation (8) 
can be obtained by a method other than the step method 
from the integral given by 


BCL, "ida 4( ВС ded ae 
Ts Jodë dy 





(9) 


whieh can be integrated numerically by taking the graphical 
slopes of В Спа Каз. 


272483—54—_1T 








Spanwise center of pressure, induced drag, and pitching 
moment.—The integration formulas of loading due to flaps 
required to obtain spanwise center of pressure and induced 
drag are not developed here. However, the integration 
formulas given in reference 1 to obtain these characteristics 
can be used with acceptable accuracy for the case of loading 
due to flap deflection. 

No rigorous procedure for determining pitching moments 
can be developed from this method. As noted in reference 
1, however, a good approximation for a wing without flaps 
can be obtained if the wing sweep is large since moments due 
to shifts of spanwise load overshadow those due to shifts of 
chordwise load. Where the effect of flaps is to be considered, 
as in the present case, it appears unlikely that the shifts in 
chordwise loads can be safely neglected. Therefore an 
estimation of pitching moments due to flaps should not only 
include the moment due to spanwise load redistribution, 
which this method will give accurately from the longitudinal 
moment of the center of load, but also should inelude an 
estimation of the moment due to chordwise redistribution. = 

Additional considerations.—Several items pertaining to 
the usage of the method should be considered. 

Ë Additive nature of loading and spanavise angle of attack: 
The linear relation between angle of attack and loading 
distribution of equation (1) states that all loadings are 
additive if the respective angle-of-attack distributions аге. 
added. 

Thus, the aerodynamic coefficients that result from the 
integration of the loading distribution are additive if they 
depend linearly on the loading distribution. Hence, Ш. 
coefficient, rolling-moment coefficient, and pitching moment 
are additive; whereas induced drag, spanwise center of 
pressure, and aerodynamic center are not. 

The additive concept is very useful for the determination 
of loading due to flap deflection. The loading due to an 
arbitrary span flap having arbitrary position on the wing 
semispan is simply found by adding or subtracting known 
loadings due to flaps at other locations on the wing. For 
example, if the spanwise ends of the flaps are at the same 
span stations, then for full-chord flaps 


C iss сага = Gi =, С outboard 


(Сайыша = Cze 一 ( Calia (1 0) 
1— 7, овеовга 

2. Differential flap angles: The effect of a differential 
between Нар angles can be taken into account by considering 
Cr, of each wing panel as having one-half the value of the 
symmetric case. The total wing lift is then the sum of the 
products of Съ 2 and the angle of deflection of each Нар. 
Although the total wing lift can be found by this procedure, 
the loading distribution cannot be found exactly by the 
products of the symmetrie 6/6 and the deflection of each Вар. 
However, the loading distribution so obtained will be suf- 
ficiently accurate for most purposes since this procedure 
neglects only the small change due to the induced effects of 
the difference in Нар defleetion on each opposite wing panel. 

3. Flap angles measured perpendicular to the hinge line: 
The relationship between a flap angle measured in a plane - 


П Паром — 
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perpendicular to the flap hinge line and that measured in 
a plane parallel to the plane of symmetry is given by 


tan è 
cos А; 





tan ô= 


(11) 


where A, is the sweep angle of the Нар hinge line, and š is the 
angle measured perpendicular to the hinge line. Equation 
(11) can be inserted into any of the equations involving 8 to 
put the equations in terms of 8. 

For constant fraction of wing-chord flaps on straight- 
tapered wings, А; 13 given by 





4[(8/4)—(c,/0)] (1-3 


tan Ar 一 tan А— 本 TEx 


(12) 


SOLUTIONS FOR STRAIGHT-TAPERED WINGS 


Charts of aerodynamic characteristies for straight-tapered 
wings can be presented in terms of parameters involving 
geometry, compressibility, and section lift-eurve slope.’ 
These charts provide a ready means of obtaining data 
direetly. ` 

Equation (4) has been used to evaluate the loading dis- 
tribution G/8,— (1/81) (cze/26) for a range of wing parameters 
and flap spans, and the results are plotted in figure 4. The 
loading presented in figure 4 gives the values due to full 
wing-chord Нарз. For other Нар configurations, the pro- 
cedure given in the loading distribution section of the general 
solution can be used. 

With the loading distribution determined from equation 
(4), the application of equation (6) provides lift coefficients 
due to flap deflection in the form ВС,,/Ка for several flap 
spans. In figure 5, these values are plotted against extent 
of flap span from the wing root outboard for a range of wing 
parameters. | 

As presented, figure 5 gives directly the lift due to full 
wing-chord inboard flaps for fiap spans measnred from the 
plane of symmetry outboard. The lift due to outboard 
flaps for flap spans measured from the wing tip inboard can 
be obtained from figure 5 by use of the relations of equation 
(10). For full wing-chord flaps located arbitrarily on the 
wing semispan, the lift can be obtained from figure 5 as in- 
dicated in the following example sketch: 


























Pts | 
Kai of the given “op 
ВС, 
“ар | 
: | Fiep span of 
— -一 一 一 He - arbitrary 
Й position 
i ara НЕ wing 
0 п | 40 


5 Throughout the figures, кал is the constant spanwise-section lift-eurve slope ór the average 
of a smali variation, Forlargespanwise variations of к that follow the function given in eque- 
tion (В5) developed in appendix В, the parameters В.А/ка, and А сап be replaced by the 

BACHA к: i 
parameters and - 4, respectively. 
(ke xd) ща 

For large spanwise variations of « that do not follow the curve of equation (В5), the simul- 
taneousequationsforthe genera! solution can be solved for arbitrary distributionsofa, Values 
of Н, can be obtained conveniently from figure 2, 
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With the full wing-chord values given above, the lift due 
to constant fraction of wing-chord and flaps of arbitrary 
spanwise chord distribution can be found through use of 
equations (7) and (9) with the da/dé values of figure 3. 
APPLICATION OF THE a,, COEFFICIENTS TO THE PREDICTION 

OF DOWNWASH FROM ARBITRARY WINGS 

In the subject report and in reference 1; it is shown how 

various forms of symmetric loading can be found quickly 


through use of simplified lifting-surface theory and values of . 


the а,, coefficients which are presented in graphical form. 
It can also be shown that if the symmetric load distribution 
is known, from any computational method or from experi- 
ment, then the same a,, coefficients can be used to find cor- 
responding values of the downwash in the center of the down- 
wash wake. It is the purpose of this section of the report 
to make available the method for computing downwash 
without attempting to explore the possibilities or limitations 
of the method. In view of the importance of being able to 
predict downwash, it is considered that inclusion of this mate- 
rial herein is justifiable although not directly related to the 
main purpose of the report. It is important to remember 
that the method can be used for any case of symmetric load- 
ing and is not restricted to that discussed in this report. 

Ав noted previously, the a,, coefficients provide a very 


simple means of finding the downwash in the center of ес. 
- downwash wake of the wing due to any symmetric spanwise 


loading distribution for any sweep angle. The downwash 
can be found at specifie points in the region between the 


wing tips from near the quarter-chord line to an infinite 


distance downstream. | 
Equation (1) was derived (reference 1) to find the down- 
wash at the effective three-quarter-chord line for which the 
distance downstream from the quarter-chord line is к(с/9). 
The general expression for the distance downstream from 

the quarter-chord line is 
г--Ву tan Ag 


where ¢ is the longitudinal distance measured in the hori- 
zontal plane from the wing-root quarter-chord point to an 
arbitrary downwash point. Substituting this general ex- 
pression in equation (2) gives | 


а, 


Н ут tan А 


(13) 


| where 


т 
Еър 


. Equation (1) (for values of the parameters ЕВ, 7,, and tan 


Ag) gives the downwash angle as 


1р 4 
(ў). = 2 0 6, 
в? n=l 


where the d,4 coefficients are obtained from figure 1 with the 


(14) 


“Н, values given by equation (13). 


It should be noticed that, because the chordwise loading 
is concentrated at the one-quarter-chord line, equation (13) 
is independent of aspect ratio and spanwise distribution of 
wing chord, while from equation (14) it can be scen that 


П 


downwash depends on these wing geometric values only 


through the loading distribution G,. 
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Equation (14) vill give the downwash angle due to arbi- 
trary symmetric loading distribution for all values of i" 
tan As. The downwash can be found for the specific span 
stations corresponding to »— 1,2, 3, 4, or n,= cos 092, 
0.707, 0.383, and 0, respectively. И compressibility effects 
are ineluded, the loading distribution must be that at the 
Mach number considered. 

For convenience, equation (14) is rewritten such that 
downwash angle and loading distribution are in terms of 
unit lift coefficient. "Then 


Low V \ 4 
(s пре га) е o). 


N Kas 


(15) 


where се Съса is the span loading distribution for unit lift. 
Both a,, and c;e/Czc,, are affected by compressibility, being 
functions of the parameters В.А/қ, and Ag. 

Equation (14) or (15) gives the downwash, in the vertical 
center of the wake, due to the integrated effect of the loading 
of а swept vortex and the trailing vortex sheet extending 
baekwards to infinity. Any rolling of the trailing vortex 
sheet into two trails is neglected. Since either of the equa- 
tions are for arbitrary loading distribution, the loading 
coefficient сап be set equal to a constant to duplicate uni- 
form loading and the downwash results obtained for the 
case of the trailing sheet rolled into two trails. 

With се/Сьсь==1, equation (15) becomes 


(7) 4 











E | 
(де Doro (eh 9 
where 
vat 
© 7 319 
‚___У 
1 ЕЗІ 


and б” (as used in the past from line-theory concepts) is the 
span between the two trails given by the integrated loading 
divided Бу the loading at the middle of the wing, or 


ыы 5 я 
Чес ерес” 


DISCUSSION 


INFLUENCE OF WING GEOMETRY ON FLAP EFFECTIVENESS FOR 
STRAIGHT-TAPERED WINGS 
To show the effect of wing geometry and compressibility 
on lift due to inboard Нар deflection, the data of figure 5 
have been cross-plotted and presented in figure 6. These 
effects are similar to those on wing lift-curve slope shown in 
reference 1. Figure 6 (а) compares the present theory with 
that of reference 3 and shows, as might be d that 


e! 


Figure 6 (b) shows the e of taper ratio on в. da: 
tiveness for swept wings. Figure 6 (c) gives the relative 














effects of the wing parameters of sweep, aspect ratio, and 
taper ratio. 

Figure 7 shows comparative effects of inboard and outboard | 
flaps on swept wings and shows the loss of flap effectiveness 
due io sweepback is less for inboard flaps. The relative 
effects indicated by figures 6 and 7 apply equally well for 
constant fraction of wing-chord flaps since the lift is only re- 
duced by a constant spanwise factor да/98. 


LIMITATIONS OF THE 2/05 CONCEPT AND COMPRESSIBILITY 
CONSIDERATIONS 


A minimum value of aspect ratio for which two-dimen- 
sional values of de/d8 сап be used with accuracy has been 
determined for use with the subject report from principles of 
lifting-surface theory and from experiment. An analysis 
based on lifting-surface theory indicates that the da/dé of a 
wing does not vary appreciably with aspect ratio until aspect 
ratio becomes less than approximately 2; for A<2, do/dó 
increases rapidly with decreasing aspect ratio until at A=0, 
da/dé=1. The discussion of reference 2 concludes that only 
& lifting-surface theory that accounts for changes in chord 
loading due to flap deflection will give accurate loading pre- 
dictions for low-aspect-ratio wings. 

Low-speed experimental values of 4о/06, from я 7- by 
10-foot wind tunnel (reference 5), are plotted in figure 8 аза 
function of aspect ratio for an unswept апа а swept-back 
wing. Figure 8 indicates the order of variation of da/dé with 
aspect ratio. 

The compressibility considerations iu the present report 
аге the same as those of reference 1, 2, or 3; namely, that 
compressibility can be taken into account by substituting 
for the geometric values of A, А, and b, the parameters ВА, 
As, and 66, respectively. 

From two-dimensional theory it can be shown that the 
parameter da/dé is independent of Mach number. In the 
present report da/dé is also considered independent of Mach 
number апа, further, is considered independent of aspect 
ratio for 4122, as discussed earlier. However, from the 
linearized compressible-flow theory the aspect ratio becomes 
an effective aspect ratio given by ВА which approaches zero 
аз M—41.0. Thus the use of two-dimensional da/dé is re- 
stricted to the Mach number range for which BA>2. 

Figure 9 shows the flap-effectiveness comparison between 
theory and experiment for a subsonic range of speed for two 
wings. The theoretical values were obtained using da/dé 
from figure 3 and assuming it independent of Mach number. _ 
The curves are terminated at the Mach number where 
BA=2. The parallelism of the theoretical and experimental 
curves of figure 9 indicates that 4о/46 can be used with the 
assumption that it remains independent of Mach number to 
values of 84 of the order of two. 

In view of the foregoing discussion, it is believed that 
two-dimensional values of da/dé can be considered independ- 
ent of M as long as ВА >2 in the application of the subject 
procedure., 


COMPARISON OF THEORY WITH EXPERIMENT 
Experimental values of fiap effectiveness are compared 
with theoretically predicted values in a correlation diagram 


given by figure 10. Included are the results for a wide range 
of plan forms. Sketches of the plan forms and flaps are 
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drawn about the points of correlation. The theory makes 
use of the curve of figure 3 giving do/d$ for a scaled-gap Нар 
over a range of deflection of +10°. Ав seen in figure 10, the 
theory is in fair agreement with experiment. 

In general, the theoretical values of figure 10 are too low 
&nd center approximately about the minus 10-percent line 
of the figure. А similar effect is shown in figures 8 and 9. 
This almost constant difference between tlie theoretical and 
experimental results is possibly due to the 4а/46 values that 
were used in obtaining the theoretical results and which were 
obtained from the —10? to 10? sealed-gap curve of figure 3. 

The theoretical flap effectiveness of the low-aspect-ratio 
wings having partial-wing-chord flaps, given in figure 10, fall 
below the minus 10-percent line. Part of this discrepancy 
is presumably due to using two-dimensional values of da/dé 
since these wings approached the lower limit of aspect ratio 
for which the da/d8 concept has been judged acceptable. 

Comparison of experimental results with theory for a Mach 
number range is shown in figure 9. This figure shows that 
compressibility effects are accurately predicted up to Mach 
numbers corresponding to values of ВА=2. 


CONCLUDING REMARKS 


The determination of symmetric loading due to Нар deflec- 
tion for arbitrary wings is shown to be obtained by the solu- 
tion of four simultancous equations. The coefficients of the 
simultaneous equations are presented in charts as functions 
of parameters that include wing geometry, compressibility, 
and section lift-curve slope. Integration formulas are given 
for the loading distribution. With these the flap effective- 
ness can be found. 

For the special case of straight-tapered wings, the loading 
distributions and values of flap effectiveness are given in 
chart form for a range of wing plan forms. These charts 
indicate that the effects of plan form on the loading distri- 
bution and lift due to flap deflection are similar to that on 
additional loading. 

On the basis of experimental comparisons, it is concluded 
that the method of the subject report can predict the flap 
effectiveness of wings at subsonic speeds with good accuracy. 


AMES AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
Мовветт Кер, CALIF., Sept. 22, 1950. 


APPENDIX A 


DEVELOPMENT OF EQUATIONS FOR DETERMINATION OF 
SYMMETRIC SPAN LOADING AND LIFT DUE TO FLAP 
DEFLECTION 


ALTERED BOUNDARY CONDITIONS FOR FLAP DEFLECTION 


In the following development the flap chord is considered 
equal in length to the wing chord; thus a deflected flap 
results in a discontinuity in the spanwise distribution of 
angle of attack. Spanwise loading for a discontinuous dis- 
tribution of angle of attack can be found through the simul- 
taneous solution method of equation (1). The number of 
equations, however, must be greater than four to obtain 
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‘accurate results and then the computation efforts become 


excessive. (Reference 1 gives the method for an arbitrary 
number of equations.) However, an alternate procedure, 
originally developed in reference 2, provides sufficient accu- 
racy with a reasonable number of equations and this proce- 
dure is the basis of the following development. 

The zero-aspect-ratio theory of reference 3 shows that for 
а given sweep angle, as aspect ratio approaches zero, span- 
wise loading becomes independent of plan form. Thus, the 
altered boundary conditions (see Development of Method) 
required for the present theory to give spanwise loading for 
given flap spans can be determined from the loading given 
by reference 3, which applies to all plan forms. These 
altered boundary conditions are obtained from equation (1) 
using the a,, values for zero aspect ratio and the loading dis- 
tribution for zero-aspect-ratio wings given by reference 3 
for given flap spans. | 

The zero-aspect-ratio values of a,, are given by figure 1 
for H,—0 or are given by the B,, coefficients of reference І. 
These a,, coefficients are tabulated as follows: 


TABLE А1 
Gin (FOR m=7 AND А--0) 







—2. 0720 0 











10.4524 
—8. 8284 5. 0568 —1.8284 
0 —2. 3888 4.3296 一 3. 1548 
—. 2928 —1.7072 " 





With these ду, coefficients and equation (1), symmetric 
spanwise loading can be found for zero-aspect-ratio wings. 
Ав а comment on the accuracy of the simplified lifting- 
surface theory for determining additional loading with 
m=7, the solution of equation (1) with A=0, a,, values gave, 
to the accuracy of three decimal places, the elliptic loading 
distribution characteristic of zero-aspect-ratio wings. 

The flap end is arbitrarily chosen for the present theory 
as the mean value of the spanwise trigonometrie coordinates 
of the points at which the boundary condition is applied and 
which bracket the discontinuity or flap spanwise end. Thus 
for m=7, three flap spans can be defined for both inboard 
and outboard flaps. 
wise point at the end of the flap, then 


1,—€o0s 6 for inboard flaps 
nr=1—cos 0 for outboard flaps 
These flap spans are given in the following table: 


| TABLE A2 
DEFINED FLAP SPANS . 


























Inboard Outboard 
pues m paren Не 
Casa I H HI IV у ҮІ 
6 Тж Би 3r Зх 5х Тк 
16 16 16 16 16 16 
| А = TEE 
т 0. 1951 0. 5556 0. 8315 0. 1685 ‚ 0.4444 | 0.8049 








Let 4, be the flap span, and 8 the span- 
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For the flap spans listed in table A2, the exact span load- 
ing distribution for zero-aspect ratio can be found from 
reference 3. With the а, values listed in table АТ and the 
exact values of Gi, б», бз, and G, determined from reference 
3, equation (1) gives the boundary condition or twist re- 
quired for the present method to produce the loading dis- 
tribution for each case listed in table A2, that is, 


5=10.4524 1.8.8284 6 € 9.2928 € = 

I H 

== 2. 07 20 S45, 6568 $a. 3888 Z da 

1 

(A1) 

lt — 1.8284 € 43296 Ф 17022 € 

6; 5; à, б 

£= —0.2242 98.1548 È sti 0000 = = 


Тһе spanwise loading distribution due to Нар deflection 
is given by (reference 3): 


Inboard Варз 











045% 
GS) sin g—4] 26 sin e--(cos е— сов 8) In] i+ 
| sin È 
2 
cos E? 
{cos p+cos 0) In (A2) 
0 一 9 
Ба "29. 
Outboard flaps 
Ез =sin = е) 
outboard à EM 


With equation (А2), the spanwise loading (1, бс, бз, and 6, 


at span stations e= 7 = and 


0.3827, amd 0 can be tabulated for each of the cases given in 
table A2 as follows: 


2 or n=0.9239, 0.7071, 


ТАВЬЕ АЗ 
SPANWISE LOADING FOR А=0 


























= Сазе 1 п ПІ ІУ ү УТ 
| E 0. 0505 0.1524 0. 2586 0.1241 0.2303 0.3322 
È 1106 3467 6432 - 0639 3604 5965 
| È 2062 -6903 8714 0525 2336 7177 
| È +4125 » 7970 + 9618 » 0382 « 2030 .5815 

















The twist distribution required for each case is obtained 
with equation (A1) and table A3, and is summarized as 
follows: 
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ТАВЬЕ А4 
TWIST DISTRIBUTION, = 
1 















































5 — 
Inboard Outboard 
Case I п п IV У УТ 
тц 
8 Е Sr SE SE dr e ad 
16 16 16 16 16 18 
я, 0.1951 0.5556 0.5315 0.1685 0.4444 | 0. 8049 
4. 4 
= —. 0170 0318 —. 0414 1.0414 9682 1.0170 
a 0285 —.0034 1.0210 —.0210 1.0084 | 9715 
+ —.0136 9041 9548 0452 . 0059 1.0136 
— = 
va | . 9884 9763 1.0404 —.0404 .0237 .0116 














With the twist distribution given by table А4, equation 
(1) can be used to solve for spanwise loading due to flaps 
for any of the six cases and any aspect ratio. 

LIFT DUE TO FLAPS 


The loading distribution resulting from flap deflection is 
distorted to such an extent that the usual quadrature 
formulas are not sufficiently accurate to integrate the loading 
for m=7. However, the following development results in а 
method which will give the lift due to flaps аз а function 
of the four loading coefficients: 


Lift coefficient is given by the integral 
1 т 
CA | am@di=A ( Gle) sin ede 
=i с 
If it is assumed the loading series 


= 
С(е)- >) а sin ие 
ыт 


then 
с," а, (Аз) 
If the loading series is expanded for e= г 2 == and : > 


апа С, Gn, Gs, and С, are obtained іп a series of а, 78, о. 
multiplying each series through by the first coefficient of the 
series and summing the G's gives 


= (0.3827 G,+0.7071 G,4-0.9239 въ 8 + 


(А4) 


The high harmonic coefficients сап be put as factors of 
the G,. With equation (А4), the lift coefficient of equation 
(АЗ) becomes 


Эд LE .3827 x (+ 


0. поте (148 ; 
0.92397 i5— Сат в) 
4 (144 1.5068 6; Ger 


т (1 Gist rr ба — 433 
"(5+ 1.5068 б, ја] 


Gis lirt da izt · + è 


ај) G+ 
1.5068 6, :十 


ar TL 6,4 
1.5068 6, 5 





(А5) 
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Н h, is defined as the coefficients of Ga, then the. lift 
coefficient is given by 


Съ А (AG, 4- ha GS t ha G, Мб) (А6) 
Reference 8 also gives loading due to flaps in series form 
from which the ratio of (aa 一 ay? 十 aa 一 ass) to 6, can be 
evaluated. It is expected that this ratio will not vary appre- 
ciably with aspect ratio. These high harmonic coefficients 
are small but are not negligible for loading due to flaps. The 
а, values from reference 3 are given by the following equa- 
tions (for odd values of ш): 
For inboard flaps 
разе 


a =a п--20--з 20) 


ш21 
46 


— sin 0 cosu, 0—cos 0 sin u;0) 
пики" — 1) 


da 


For outboard flaps (for the same 0) ay) 


pic 
(41) outdosra = б— (а 1)inbosrd 


ul 


(ал) outboard = 一 (а,) inboard 


With equations (A5) and (A7) and the values of table A3, 
the А, values for each case can be tabulated: 


TABLE A5 
LOADING INTEGRATION FACTORS, ha 





























| as У п ш | | ~ | vol vi | 
z h | оз 0, 2994 0.2999 | 0.3020 0. ша. rl; 3008 
“ЕТ . 5541 . 5544 3549 5608 з ЕЕЕ 

мо. 27248 „7250 27252 21339 sums |) .ю% | 

м .3922 .3921 „3922 | 4080 | .3950 .3990 | 








INTERPOLATION FUNCTION FOR SPANWISE LOADING DISTRIBUTION 


The method of this report gives spanwise loading due to 
Нар deflection at four span stations. Since the complete 
loading distribution is not sufficiently defined by these values, 
the values of loading at other span stations cannot be inter- 
polated accurately by direct use of the interpolation equation 
and table of reference 1. 

The direct use of the interpolation table of reference 1 can 
be obtained by converting the present loading due to flap 
deflection to a loading approximating additional loading dis- 
tribution. The conversion factors are obtained from values 
of the zero- -aspect-ratio theory and given by the ratio of 
equation (А2) to sin о, for each of the four giv en span 
stations. Define 

__(G/8)4=0 
sin e 


then, with table АЗ and equation (АЗ), values of А, have 
been computed and are tabulated for each of the flap span 
cases. 


(A8) 
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TABLE Аб 
Iuboard Outboard ” 
Case I n у І v 
* т $ | s 
n | 0.195 0.556 0,169 0.54 
Я 
E N БЕСТЕ М 
1 0.924 | 0.1318 0.3081 0.3243 | 0.6019 
201600 1564 4903 .9041 5097 
= > RENNES 
3 - 383 2232 ‚Тїї 5685 . 2529 
4 0 #25 | „0 0382 . 2030 

















where, for 7;=1, R,—1.- 
The interpolation table of reference 1 gives 


TABLE A7 
INTERPOLATION FACTORS екі 












































The interpolated loading at span station Ё is determined 
from the summation 
4 
то экы ein 
i Е Бела 
where oqtiations (A2) and (A8) аи values of Ry which | 
are tabulated as follows: | 
TABLE AS. 


- (А9) 








ааа ay, ы» aisi 
ГА 0.9808 0.8315 | 0.5558 | 0.198 
` | 4 
NE 1 3 б nt 
" 2 3 3 vg 
N | ИРИНА 
1 | 0.8657 0.4904 |.-0.0976 | “020229 Е 
Pale Маа 
2 | „5130 | —; 0746., 
PRA МЕНЕ = 
3 | тп | .4187 
























































„2 Р, 2 
` Case 
Ер І п ш IV Y 
7% 
— — 
š ' 0. 981 0.1246 0.3826 0. 6090 0.3911 1 0.6176 0.8754 
m === ; Е 
3 а ^ 
3 E 831 . 1421 ‚4304 ‚8090 .1910 .5695 . 8580 
8 М + 
2 = 558 . 1799 . 6253 . 9314 + 0685 . 3746 .8201 
“ al 
1 .195 . 3324 . 7863 . .9556 . -0442 .2136 .6076 








With the R, values of table А8 and the loading ratios & 
k 


determined from equation (A9), с loading G can be found 


at span stations k= - 2. aut 2 (1—0.981, 0.831, 0.556. 


2 2 
and 0.195, Жете 


APPENDIX В 
AIDS TO THE METHOD 
LINEAR ASYMPTOTE OF ayn 
For large values of H,, the а,, functions become linearly 
proportional to Я, (see, e. g., fig. 1). Since thislinear charac- 
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teristic appears at relatively low values of H, the simple 
linear relation between a,, and ZT, is quite useful. 

The L,,* function of appendix A of reference 1 is multiplied 
by b/c, and the product is linearized with respect to b/c,. 
The result gives 


b 
+): == sol = nns 
2 





























— 7 25. 
jj Po Эт tan À E 
= Iob Элуу Gs X 
Aa S Sm ATT астар 
[+ = ола) —1 (В1 
\ ыш ал for 7< pen 
27 tan A xi 
= ( 二 一 一 tan ^) (2) 
с е, n 
1 1 А = 
an Тэг 1 5 for 7= 7 
PM ) 2) 2 с 
----9 | У па кан 
2 безь ўтва А (=. + (1--віп А). 
for п=0 


With the values of equation (B1) substituted into equation 
(A37) of reference 1, the values of a,, for arbitrary sweep 
angle and for high values of H, are obtained. For m=7; the 
following equation (B2) gives values for a,, as. 


ЭЕ (2:313 .913 


+1.452 tan а) H+ ава: 751 sin A— 
cos А ' 





/1--0.0016 tan? = 
У — 
0.003 tan = 5+0. oss ( mm 


TL 77 2A 
сова (= - 0.0177 tan? А De 
































tan À 
0.014 (is tan? ^=) 
' [ап А 
и (В 3.798 {ап А) Н,-3.100--1.610 sin A+ 
г ha 2 人 _ | 
0.008 tan s—0.048 («рые Se 1). 
2 tan А 
+0.0177 tan? A—1 ve) 
мая (95 tan А )- 
+0.1716 tan? A—1 
2 
OAs (93 бап А ) 
2-0.25 Я 
ала 29 E2416 tan A) H,—0.014—1.727 sin A— 





/ 





/{+0.0016 tan? zr 


0.032 tan 340. 037 (18 E + 


PTS Шен tan? E 








tan À 





1140.1716 гс, и 


E 
LOS ( tan А 








_ (0.084 


(22820. 408 tan a) H,—2.140—0.293sin A— 


”Ca 一 


0.067 
— же 3E 
Se ( соз А. ! os 


0.079 ава + Me 





0.060tan_0.017(* 








i1+0.0016tan*A—1 











0.048(* 
0. oss(^ 


0.004tand0.003(* 


тап А 





— A1) 


tan А 


tan А 


кыр s) 








влл1в(У 
oo1s (A12 





voti 


/ 
0.063(* 


0.088(* 


0.275 
cos A 


0.042t3—0.107 (XL EGIT Ian 


10.0294 tan?A— D. 
tan А 





tan A 


+0.0886 — 1) 








tan А. 


1--0.0886їап?А—1 


1+0.0294tan?A— 1)- 


tan A ) 








0.048(* 


tan А 


414-0. аи) 








0.063(* 





y1+0.0886tan? =) 


tan А 





0.022 (21 


^12-0.0294tan*A—1 
J+ 


{ап А. 








o.116(* 


0.041 


—————0.068їапА ЈЕ 


cosà 


0.007 tan*—0.1 16 


ni 
0.125 (° 


0.2134 (2 


1+0.0886tan? ы) 
tan А 


—0.043+0.637sinA— 


)- 


1--0.017 Tian 4-1) 
{ап А. 


tan У 4.640— 0.491 sin A+ 


414-0.0177 tan? х1) 
tan А 


— 0.986 tana) Н.-2.9264-0.8218щА-- 


14 шкы 
tan А 


7 tan а) H,+0.188—0.551 sin A— 


1--0. ТОТИ) 
tam А 


0.197 4- 0.207 sinA— 





tan А 


/1+0.1993 tan?A—1 
tan A 





) 


(4 +0.1716 tan?A— D- 
2 tan А 


1+0. FOSSA anti 1). 


i 
| 
ü 
| 


255 


(В2) 





256 





по ШЫЛ 0.281 tan л) 太一 1.850 一 0132sinA 十 








41--0.1716tan?A =1).. 


A È 
0.020tan3+0.302(* tan A 








T (А --0.0886 tan*A— b. 
: {ап А 








, La 
of 14-0.1993 tan 1) 
tan А 





(0:522 + 0.282 шпА)Н,4-4.272 —0.388sinA— 
cos A 








— 1) 


= A mf x 
0.078tan3—0.197( IDA 


/1+0.0886 ian Ai) 
š ` 
0.302 (| ui E 











і 
0.088 ( ni 


0.156 | 
a= (50.428 tan A) H,—1.666--04895in A— 





| / 7 YA nc 
0.145tan 人 十 0.088 41--0.1716tan?A 1)- 


tan Á (B2) 
sa NEL A ў 
IIS e 
tan А 
0.041(V1+0.1993tan'A_1 
tana 


250.032 (tan A) И, 0.224 
1 ; 
dag=—0.097(-2+tan ^)H.-0.354 (1--віп А) 
1 
аз=0.371 (acit) Н;—3.155 


1 | 
240.090 ( tan) He+3.313—0.687 Sin А 


К. 


LINEAR SPANWISE DISTRIBUTION ОЕ (кс), /(кс) е 


Equation (3) сап be manipulated into a function of various 
parameters involving distributions of wing chord and section 
lift-curve slope. One procedure is as follows: With the con- 
dition that the product of section lift-curve slope and wing 
chord varies linearly spanwise, then 

(1-8 221 


LS: 
2[1 —3,(1 — к.Мк)] 


where A, is the aspect ratio based on a wing chord equal 
to (xc). Equation (B3) is reduced to terms of two param- 


Cose ere rev) [77 38 


апа equation (3) becomes 


Н,=4,(8А) 
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eters, A, and (к/к)Х. Expressions of A, in terms of aspect 
ratio for straight-tapered wings, and the distribution of sec- 
tion lift-curve slope can be found. 

For straight-tapered wings 





2b 
А-а») айы 
and since (ке) is linear 
4 2b i Ее к 
548 krCr [1 + (к/к) 
Шеп 
ne к. K À 5:5 m 


and equation (B3) becomes 


_3 Гвай+ж) 14 (k/x)A 
Нина, | М), (1—7, 1 —(кјк ЈА] ) 





“ва 


The distribution of к for straight-tapered wings is given by 


CIA 1— n» [1—(кјк)Х] » · (B5) 


Су * 1—3,0 —) 








EU 
(ке кА) 
апа (к;/к,)А. Solutions for spanwise loading in terms of these 
two parameters and Ag are valid for the distributions of sec- 
tion lift-curve slope given by equation (B5). Equation (B5) 
indicates that at X —1, к, is a linear function and at \=0, к, із 
a constant. For values of X between 0 and 1, x, is a curve in 
the region between the linear function and a constant. 
Equation (B4) is represented in graphical form by figure2 but 
IL 
ВАСЕ Ки кол) 


Equation (B4) is in terms of the two parameters | 


with the ordinate replaced by the parameter 


and the abscissa by (к/к. 

In summary, for straight-tapered wings having & linear 
distribution (xe), the loading and associated aerodynamic 
characteristies can be presented in terms of the parameters 

BAHA) 

P a Fk 


з and (кұқы 
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TABLE I—SYMMETRIC INFLUENCE COEFFICIENTS, а,» 
BEYOND THE SCOPE OF FIGURE 1 
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FIeuRE 1.—Influence coefficients, бал, for symmetric spanwise loading plotted as a function of the wing-chord-distribution parameter, 7 Г, for values of the sweep parameter, Ag, degrees, 


(8) veil, пе 0.061. 
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«FIGURE 1.— Continued. 
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FIGURE 1.— Continued. 
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FIGURE I.—Continued, (g) »=2, п=3, 41=0.234. 
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FIGURE 1.—Continued. (k) v=3, п=3, d:=0.381. 
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Figure 1.—Continued, (0) уш4, n=3, di=0.320, 
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Fieure 2.— Variation of the geometric parameter, x with taper ratio, А, for straight-tapered wings, 
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FIGURE 3 — Variation of lift-effectiveness parameter with fiap-chord ratio, z. Average trailing-edge angle about 109; М=0.2. (Curves from reference +.) 
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FIGURE 4.—Spanwise load distribution, E per radian, due {о symmetric flap deflection for straight-tapered wings. Z= 1. 
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FIGURE 4.— Continued. 
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FIGURE 4— Continued, 
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FicuRE 4.—Concluded. 
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FIGURE 5.—Flap lift parameter, per radian, for extent of flap span from the wing root outboard, ч =1. 
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FIGURE 9.—Comparison of theoretical and experimental flap-effectiveness увла ог with Mach number. Reynolds rumber=2X10% at М=0,2; 4.8X10* at Af=0.6. 
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FIGURE 10.—Correlation of theoretical and low-speed experimental fiap-effeciiveness due to symmetrically deflected flaps for various plan forms. Reynolds number from 1to3XI0 Af<0.2. 


